Fermi-Liquid Theory for Anisotropic Superconductors 
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This article develops a Fermi-liquid theory for superconductors with anisotropic Fermi surfaces, 
Fermi-liquid interactions, and energy gaps. For d-wave superconductors, the Fermi-liquid interac- 
tion effects are found to be classifiable into strong and negligible renormalizaton effects, for sym- 
metric and antisymmetric combinations of the energies of k f and — fc J, quasiparticles, respectively. 
Furthermore, the leading clean-limit temperature-dependent correction to the superfluid density in 
a d-wave superconductor is found to be renormalized by a Fermi velocity (or mass) renormalization 
effect. The question is raised of whether or not the penetration depth in the high temperature 
superconductor Yl^CusOe+^a can be accounted for with physically acceptable parameters within 
the framework of a quasiparticle model. Fermi-liquid corrections to the spin susceptibility and to 
the zero-energy magnetic-field-induced density of states are also evaluated. 

PACS numbers: 74.20.-z, 74.25.Jb 



I. INTRODUCTION 

There is now considerable experimental evidence that 
the cuprate high T c superconductors exhibit the sim- 
ple power law temperature dependences predicted by the 
quasiparticle picture for their thermodynamic and trans- 
port properties at temperatures well below T c . For ex- 
ample, penetration depth measurements find that the su- 
perfluid density exhibits a low-temperature clean-limit 
linear-in-T temperature dependence^, in agreement with 
theoryB. The NMR relaxation rate exhibits the ex- 
pected T 3 temperature dependence^. The predicted ef- 
fect of impraities in giving rise ta a universal thermal 
conductivity^] has been confirmedta. The clean-limit spa- 
rine heat varying as T 2 appears to have been observedEra. 
Even the electrical transport relaxation rate observed in 
microwave conductivity experimental, which had resisted 
explanation for some time, hap-Jiow been explained in 
terms of a quasiparticle pictured. 

Whether or not the Fermi-liquid parameterization of 
the coefficients of the above power law temperature de- 
pendences is quantitatively accurate is at present an 
open question— A recent study correlating these differ- 
ent coefficients^ concludes that the quasiparticle model 
may be successful here also provided a Fermi-liquid inter- 
action factor multiplying the superfluid density is treated 
as an adjustable parameter. Some remarks at the end of 
this article address the question of whether or not the ex- 
perimentally determined value of this adjustable parame- 
ter has a physically reasonable value. The answer to this 
question should help to assess the validity of the quasi- 
particle picture of the low temperature properties of high 
T c superconductors. Recent debate on correctness of the 
quasiparticle picture_iS|-a]so occurring in connection with 
ARPES experimentsHLia, and in connection with the role 
of phase fluctuations of the complex order parameter in 
the determination-of the temperature dependence of the 
superfluid densityllJ. 



A principal goal of this article is to develop a 
Fermi-liquid theory applicable to superconductors with 
anisotropic Fermi surfaces, Fermi- liquid interactions, and 
energy gaps, and to apply it to the c£-wave superconduc- 
tors at temperatures well below the critical temperature 
T c . This theory will then be used to develop formulae, 
including Fermi-liquid corrections, for the London pen- 
etration depth, the spin susceptibility and the magnetic 
field contribution to the density of states and specific heat 
(observable in the mixed state) . The approach used here 
is a relatively elementary phenomenological ope, similar 
in spirit to that of Landau's original articlelij and its 
extension_io superconductors by Betbeder-Matibet and 
NozieresEEl, and it is hoped that it will be of pedagogi- 
cal interest. An important aspect of the extension of the 
Landau theory of Fermi liquids to superconductors is the 
introduction of an appropriate dependence of the energy 
functional on the. superfluid momentum. The classic ar- 
ticles of Larkintll and of LeggettEa used a more formal 
correlation function approach than is used here. Other 
studies of Fermi-liquid interac tions i n unconventional su- 
perconductivity include Refs. 
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It should be emphasized that Fermi-liquid theory is 
known not to give an adequate description of the normal- 
state properties of the high T c superconductors. Thus, 
in the application of the results of this article to high- 
temperature superconductors, it is only the properties of 
the superconducting state at temperatures well below the 
critical temperature T c that will be considered as being 
possibly explicable in terms of Fermi-liquid theory. 

The potential importance of Fermi-liquid interactions 
in renormalizing the superfluid density in the case of 
d-wave superconductors has been emphasized in Refs. 
^l| and ^2]. The calculations of Fermi liquid proper- 
ties in these articles were carried out within the frame- 
work of models with isotropic Fermi surfaces and Fermi- 
liquid interactions, but with d-wave energy gaps, and 
they demonstrated the existence of a factor renormal- 
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izing the superfluid density that was given in terms of 
the isotropic Fermi- liquid interaction parameter F{. In a 
preliminary account of the present worked a similar factor 
renormalizing the superfluid density was found, but ex- 
pressed as a ratio of two different types of velocities (i.e. 
roughly speaking as an effective mass ratio). By using 
the Landau effective mass relationE-3 known to be valid 
for Galilean invariant systems, these two results could be 
seen to be equivalent, but the relation between the two 
results for the case of systems with an underlying crys- 
tal lattice was far from clear, and has been the subject 
of some discussion and controversy. Similarly, the ex- 
pression used to calculate the current (analogous to Eq. 
p4| below) in Ref. ^] nas been controversial. For these 
reasons, the basic ideas of the approach to the Fermi- 
liquid theory of anisotropic superconductors used in this 
article are outlined in some detail in Sections II and III 
below. The formula developed below for the penetra- 
tion depth is shown the be in agreement with the re- 
sults of LarkinED and of Leggettll3 for s-wave supercon- 
ductors with isotropic Fermi surfaces and Fermi-liquid 
interactions, and with those of Refs. ^l] and |2^ for d- 
wave superconductors with isotropic Fermi surfaces and 
Fermi-liquid interactions. The result below for the low- 
temperature penetration depth for d-wave superconduc- 
tors is also valid for the case of anisotropic Fermi sur- 
faces and Fermi-liquid interactions. A formula applica- 
ble to this case has also been developed independently 
and stated in Ref. |24|, which finds a result expressed in 
terms of parameters that are different from those in our 
formula; there is thus the possibility that these two re- 
sults are not equivalent, as a result of the absence of a 
generalization of the Landau effective mass relation to 
the case of an anisotropic Fermi surface. 

The renormalization of generalized external fields by 
Fermi-liquid interactions in a d-wave superconductor at 
low temperatures has an interesting symmetry prop- 
erty. This manifests itself when the quasiparticle ener- 
gies are separated into parts that are (Symmetric and 
Antisymmetric combinations of the energies of the +k f 
and —k I states. (The calligraphic letters S and A are 
used here to emphasize the difference with the more usual 
definition of the symmetric and antisymmetric combina- 
tions with respect to +k f and +k J, states common in 
normal state analyses, e.g. see Eq. 1.32 of Ref. In 
the presence of Fermi-liquid interactions, the (Symmetric 
and Antisymmetric corrections to the quasiparticle en- 
ergies obey integral equations that are independent of 
each other, and they are renormalized differently. This 
leads to the fact that the (Symmetric external fields ex- 
hibit strong Fermi-liquid renormalization effects to which 
quasiparticles from the entire Fermi surface contribute 
(cf. Ref. ^l|), while the the Antisymmetric external fields 
exhibit relatively weak temperature-dependent renormal- 
izations that arise from the nodal quasiparticles only, and 
that can often be neglected. 

Temperature gives a (Symmetric correction to the 
quasiparticle energy because +k f and — k j states are 



affected in the same way by temperature. A super- 
fluid flow generates an Antisymmetric correction since 
the components of +k and — k along the superfluid ve- 
locity have opposite signs. Also the Zeeman interaction 
generates an Antisymmetric correction because the spin 
t and spin j, contributions to the energy have opposite 
signs. Thus the superfluid density and the magnetic sus- 
ceptibility are negligibly renormalized by Fermi-liquid in- 
teractions, while the effects of temperature (although rel- 
atively small) are strongly renormalized by Fermi-liquid 
interactions. 



II. HAMILTONIAN AND CURRENT DENSITY 

The Hamiltonian describing a system of electrons in 
a periodic lattice potential V per (r), a magnetic field de- 
scribed in terms of a vector potential A(r), and interac- 
tion through an electron-electron interaction 7ii n t is 

H = J <Pr^ a {Y)H^a{r)+H mt (1) 

where 

H = ^- (~ihV + p.) + V per (r) - fi (2) 
Am 

and 

Ps = W(0/2) - -A. (3) 
c 

Note that a term — \xN has been included in H. Also, 
the current density operator is given by 

J(r) = ±- f [(-iftV + Ps)* CT (r)] f *„(r) + h.c. j (4) 
2m I J 

where h.c. indicates the Hermitian conjugate of the pre- 
ceding term. The above four equations are probably more 
commonly written as above but with 9 = 0. The substi- 
tution ^o-(r) - * ^ a {r)e l6 / 2 in the more common form of 
these equations yields the equations as written. In the su- 
perconducting state, p s , which is a gauge- invariant quan- 
tity, will be interpreted as the superfluid momentum. In 
calculations of the penetration depth in the London limit 
below, p s can be assumed to be uniform. 

The thermodynamic potential f2 and the equilibrium 
density operator p are given by 

n = -* B nn(tre-^) J P=— W - (5) 

It is easily seen that the average value of the current 
density operator, defined by J(r) = trpJ(r), can be cal- 
culated by taking the functional derivative of the thermo- 
dynamic potential £1 with respect to the vector potential 
A(r), i.e. 

J(r) = -cSft/6A(r). (6) 
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If the current density is homogeneous, which will be the 
case of direct interest below, it can be calculated by using 
the formula 



_ 1 f , 3 ~, , e on 
J = -y^ r trpJ(r) = - — 



(7) 



The energy eigenvalues of the Hamiltonian Hq for 
p s = 0, which describes a single electron in a periodic 
lattice potential, form energy bands with the allow states 
characterized by a wave vector k lying in the first Bril- 
louin zone. Subsequent discussion will consider only a 
single energy band, whose energy spectrum is described 
by the function e b (Kk). The general problem of the en- 
ergy spectrum of an electron in an external magnetic field 
is quite complex, but for sufficiently weak magnetic fields, 
i.e. a sufficiently slowly varying vector potential, a quasi- 
classical approximation can be developed which gives the 
energy spectrum of a siagle band of electrons in the form 
e|(p s ) = e b (hk + p s )c3. The single-band Hamiltonian 
envisaged in this article thus has the form 

ka 

7777 E Vk 1 k 2 k 3 k 4 c\ lcri c\ 2rT2 Ck3a2 C k4a 1 ■ (8) 



27, 

kxk^k^k^u ier2 

The interaction term has been constructed so as to be 
isotropic in spin space (i.e. no spin-orbital interactions 
are considered), and momentum conservation to within 
a reciprocal lattice vector is implicitly assumed. 



III. FERMI-LIQUID LIQUID THEORY FOR 
ANISOTROPIC SUPERCONDUCTORS 

As an introduction to Fermi-liquid theory for super- 
conductors, a Hartree-Fock-BCS approximation is de- 
veloped. This gives an approximation for the parame- 
ters occurring in the Fermi-liquid theory which has some 
interest, and also serves to give an indication of the 
structure expected for a full Fermi-liquid model. The 
Hartree-Fock-BCS approximation follows from averaging 
the Hamiltonian of Eq. ^ over a BCS type of state, i.e. 
by assuming a state in which the excitations are inde- 
pendent Bogoliubov quasiparticles. This gives, for the 
average energy (minus /iN) 



kak' 



(9) 



where 

fkk> =Vkk>k>k~ Saa'Vkk'kk', 9kk' =Vkkk'k' ( 10 ) 

with k = — fc and 



n k <T = {c k<T Cka), W k = (Cfctc-fcj). 



(11) 



In Fermi-liquid theory, the energy is written in terms 
of the ground state energy, plus a contribution due to 
excitations from the ground state. Thus, in Eq. ^ put 
Srika, where n.2 is the ground state value of 



nko 



n 



o 



rika in the superconductor. This gives the average energy 
in the form 

e = E Q (p s ) +y^£fc -(p s )fa fcff + 

^7 E fkk'^i ka 8n k , al + ^Y^9kk'W;W k ,. (12) 



kak' * 



kk' 



Here 



£ fc(Ps) — £k 



(13) 



In the Hartree-Fock-BCS approximation the energy of 
a quasiparticle Sk in the absence of other excited quasi- 
particles, and what is called here the electrical transport 
velocity Vj., are given by 



k' a' 



kv n °k' 



l de b k 



nit- (14) 



Notice from Eq. |l4| that in the Hartree-Fock-BCS approx- 
imation the quasiparticle energy, and hence the quasi- 
particle velocity v,t = h~ 1 d£k/d]<., are changed from the 
bare band energy and band velocity by a term due to the 
electron-electron interaction, but that the transport ve- 
locity is unchanged by the electron-electron interaction. 

Consider a band kinetic energy of the form e b k — 
p 2 /2rrib (in this case the effect of the periodic potential 
is represented by replacing the free electron mass m by 
a band mass nib) so that the basic one-electron states 
are the plane waves exp(i\i ■ r). Also assume an electron- 
electron interaction v(ri — r^). Then the Hartree-Fock 
approximation to f%%, will have the form f aa (k — k'). 
In this case the equation for obtained from the first 
of Eqs. can be transformed into the relation 



Vfe = v fe-77E/ CTCT '( k - k ' 



de k ' 



Vfe- 



(15) 



The model just described has the property of Galilean 
invariance, although for particles of mass m\, rather than 
the free electron mass, and Eq. [l5]is the Landau effective 
mass relation that follows from this propertyc-3. For wave 
vectors k on the Fermi surface, v£ = fik/nib, while the 
relation v,t — hk/m* defines the quasiparticle effective 
mass m* . In the general case of an anisotropic Fermi 
surface, where the quasiparticle interaction function f%%, 
depends separately on k and k', no relation analogous to 
Eq. |l5| can be derived. By using v& = h~ 1 de k /d]s. and 
the Hartree-Fock result Eq. |lj, an integral equation is 
found that presumably can be solved to find in terms 
of and the Fermi- liquid interaction parameters if these 
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are known at all points within the Fermi surface, but it 
is not clear that this would be a useful approach for a 
more general Fermi-liquid theory. 

For models that are isotropic and Galilean invariant it 
is a known exact result (e.g. not restricted to the Hartree- 
Fock approximation just discussed) that the transport 
velocity is not renormalized by electron-electron inter- 
actions, and hence that the transport velocity is Kk/m 
where m is the free electron mass when there is no pe- 
riodic potential, or the so-called crystalline mass mb in 
a model in which the effect of the periodic potential is 
represented by replacing the free electron mass by the 
crystalline mass in the kinetic energyEJl Thus the elec- 
trical transport velocity is quite a different quantity from 
the quasiparticle velocity. 

If the electron-electron interactions are sufficiently 
strong, the Hartree-Fock-BCS approximation is not ade- 
quate. Nevertheless, the energy E of the system can be 
expected to have the same form as in Eqs. [l2] and |l5| , 
and with the parameters , gtk', £fc and v| subject to 
the constraints of symmetry but otherwise to be treated 
as unknown parameters, some of which might be deter- 
mined from measured experimental properties. Also, the 
ground state energy, which must be an even function of 
p s , is parameterized in the form 



1. 



E (p s ) = E (0) + ^Vy2S a pp aa p a p 



ap 



(16) 



An assumption implicit in this approach to the Fermi- 
liquid theory of superconductors is that the normal-state 
quasiparticle properties are not changed much in the su- 
perconducting stateta, given that the maximum gap is 
much smaller than the Fermi energy. This assumption 
is fulfilled more strongly for conventional superconduc- 
tors than for the high T c superconductors which will be 
the principal interest below. However, in the investiga- 
tion below of the properties of high superconductors at 
temperatures much lower that the maximum gap, the 
gap will be approximately constant, and this assumption 
should not give difficulties. ,— , 

The distribution function can be writtenES as the two 
by two Hermitian matrix 



nfcf 



Wt 



W k 1 - n fe j. 



(17) 



so that the variation of the energy (Eq. [12|) with respect 
to variations in the distribution function is 



5E = y^Jr(e k 8n k ). 



(18) 



Here the matrix representing the quasiparticle energy in 
the presence of other quasiparticles is given by 



£/ct(Ps) 



A fc -£fei(p s ) 
where i kl r(p s ) = s klT (p s ) + Se ka , in which 



(19) 



Seta = +i fkk' Sn k'cr>; 



also, 



fcV 



k 



(20) 



(21) 



The entropy of the Fermi liquid is 

S = -fee ^ tr[(n k lnh k + (1 - n k )ln(l - h k )] (22) 
k 

and the thermodynamic potential is = E — TS (recall 
that E includes the term —fiN). Minimizing the ther- 
modynamic potential with respect to variations in the 
distribution function yields the equilibrium distribution 
function, which is 



h k = [exp(e k /k B T) + 1}' 



(23) 



Now consider the derivation of a general expression for 
the current density in a superconductor by using Eq. [t]. 
It should be noted the the thermodynamic potential de- 
pends on p s implicitly through the dependence of the 
distribution function on p s , and explicitly through the 
dependence of Eo(p s ) and e k (p s ) in Eq. [l2] on p s . Be- 
cause the condition for equilibrium is that the derivative 
of f2 with respect to the distribution function is zero, the 
implicit dependence of the distribution function on p s 
can be ignored in applying Eq. f^. The current density 
from Eq. is thus 



P ka 



^ vl a 5n ko 



(24) 



The first term is a temperature- independent contribution 
due to the flow of the condensate, and the second is the 
contribution of the excited quasiparticles. 



IV. RENORMALIZATION OF EXTERNAL 
FIELDS - BASIC EQUATIONS 

When an external field is applied to a Fermi liquid, it 
excites quasiparticles, and these quasiparticles renormal- 
ize the contribution of the external field to the quasipar- 
ticle energy. This is the effect that must be evaluated in 
order to evaluate such effects as the screening of an ex- 
ternal magnetic field by the superconductor (and hence 
the penetration depth), and the spin susceptibility. 

The quasiparticle Hamiltonian describing the excita- 
tions of the superconducting state has the following form: 



' k r 



-kll 



Ck + Afe A fc 
A fe -C k + X k 



('k\ 

-fel 



(25) 
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where the two by two matrix is the matrix i k of Eq. 
|l9| generalized to include arbitrary external fields, and 
reorganized into (Symmetric and .Antisymmetric compo- 
nents. Here ( k = e fe + 8ef + /if, A fe = Sei + hi, and 
Afc is the momentum-dependent gap function appropri- 
ate for d-wave symmetry, which is taken to be real. An 
important step in the analysis, as described qualitatively 
above, is the separation of the excitation energies into 
Symmetric and .Antisymmetric parts defined by 



f 



1 



fc fc = o t fcfc T " Ss -kl] > Set = - [ fc fcT + fc 



■Hi 



(26) 



where interaction parameters-, ff^) (in contrast to the 
more conventionality definedcJ parameters f kk ,) are de- 
fined by 



/(i) I f uu _i_ fcra \ rs.a 1 / raa i f au \ I ot\\ 

Jkk' ~ 2^ kk ' Jk-k'Jt Jkk' ~ 2\Jkk' ^ Jkk')- l du 7 



In applying these equations to a copper-oxide plane of a 
high T c superconductor, the volume V should be replaced 
by the area L 2 . Also, /(e) = [exp(e/k B T) + l]" 1 . 

The gap equation can be found from Eqs. O, |2l] and 
I27L to be 



The quantities hi and hf in H represent generalized 
external fields. For example, the case of an external mag- 
netic field acting on the orbital motion of the electrons 
corresponds to the external field h£ = \\ ■ p s , h k = (cf . 
Eq. [13]). The case of an external magnetic field H act- 
ing on the spin degrees of freedom is described by taking 
hi — ^bH, hf = 0. In both of these cases, the magnetic 
field acts only on the .Antisymmetric mode, and has no 
effect on the (Symmetric mode of excitation. 

In addition to causing changes in the energy of a quasi- 
particle (as in Eq. , excited quasiparticles can give rise 
to changes in the gap functionEj. There are however no 
changes that are linear in the superfluid momentum, and 
this effect will therefore be neglected. (This follows from 
Eqs. |7| and |l[ See also Ref.|o|) 

The diagonalization of the Hamiltonian of Eq. gives 



ft = X! E k,all t(T lk,G, Ek,„ 



Eh.rr = E, 



IT A' 



hik 



kcr 



(27) 



v H 9kk ' 



2V 



E L 



■[1 - /(Em) -f(E-v (31) 



V. PENETRATION DEPTH FOR ISOTROPIC 
S-WAVE SUPERCONDUCTORS 

The formalism just developed can be tested by deriving 
a formula for the penetration depth in the London limit 
for the case of a not necessarily translationally invariant 
model having an isotropic Fermi surface, Fermi liquid 
interactions and energy gap, and comparing the results 
with those obtained in the classic articles by LarkinllZl 
and by Leggettlia. To do this, consider Eq. [2{j with the 
external field h k = v£ • p s describing the interaction of 
the orbital motion of the quasiparticle with a magnetic 
field. Because hi is odd in k, the parameter ff£) can 
be transformed into f kk ,. Also, for an isotropic Fermi 
surface 



A k , and the j' k 



are opera- 



where a = ±1, E k = y/Q 
tors creating Bogoliubov quasiparticles. 

Later, the energy E^ = y/£| + A 2 describing the 
quasiparticle spectrum in the absence of other excited 
quasiparticles is also used. For a d-wave superconductor, 
the quasiparticle energy can be parameterizeda in the 
neighborhood of the Fermi-surface nodal points (see Fig. 

|l|) as E^ = yj (pivf) 2 + {P2V2) 2 , where pi and P2 are 
components of the momentum relative to the nodal point 
in directions perpendicular and parallel to the Fermi line. 
At low temperatures, only quasiparticles close to these 
four points can be thermally excited. 

Using Eq. ^ in Eq. |2^ and keeping only terms up to 
linear order in the fe's and h's yields the integral equa- 
tions 



Set 



and 



V" t(+) 



1 kk' 



£fc' 



(Oh 



A 2 



h S k> ) 



Sei 



^f(-) d f 

Jkk 



V ^ 

k' 



dE 



(0) 



(5e 



(28) 



(29) 



e k = h(k - k F )v Fl f k ' k a , = V f e ' a P e [cos(k ■ k'/k 2 F )j. 



=0 



(32) 

The solution to the integral equation |2^ now gives the 
renormalized value of the external field h k — vt ■ p s as 



hi 



6ei = vi.p s /[l + lF{f(T)} 



where 



of 



Ko)^Va4 0) ' 



(33) 



(34) 



and ^(0Uf= k 2 F /(2Tr 2 hv F ). The quantity f(T) can be 
regardedEj as an "effective density of single-particle lev- 
els" near the Fermi surface relative to the normal-state 
value u(0). Note that f(T c ) = 1 and /(0) = 0. The con- 
tributions 6e k to the quasiparticle energy will be assumed 
to be negligible because of thepresence of the factor e k 
on the right hand side of Eq. g8| 

A more restricted isotropic model is now considered 
in which, following Betbeder-Matibet and NoziereslIS, 
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the band energy is assumed to have the form e\ = 
h 2 k 2 /(2rrib), i.e. the effect of the underlying periodic 
potential is simulated by introducing a crystalline mass 
771/, in the place of the free electron mass. The impor- 
tant thing here is that the kinetic energy is quadratic in 
momentum, which is a necessary condition for Galilean 
invariance. Secondly, it should be noted that the as- 
sumption of isotropy (rotational invariance with respect 
to arbitrarily placed rotation axes) automatically implies 
invariance with respect to an arbitrary translation, since 
a product of two rotations about different but parallel 
axes through equal but opposite angles is equivalent to 
a translation. Thus isotropy implies translational invari- 
ance, and these symmetries, coupled with a kinetic en- 
ergy quadratic in momentum imply Galilean invariance. 
Hence the Landau effective mass relationt3 correspond- 
ing to non- interacting particles of mass vn^ applies, which 
yields v F = vp(l + F* /3). This relation will be used to 
eliminate v F from future formulae. 

Now using Eq. |3^ together with Eq. |24| yields, for the 
inverse square London penetration depth, 

J_ 4e^ Fl T _ Ff f(T) 1 
A| 3tt?7 C 2 [1+ 3 j [ { + 3 J l + iFf/(T). ' 

(35) 

It is the second term in Eq. |24| (i.e. the quasiparticle con- 
tribution to the current) that gives the second term in the 
square brackets in the above equation (which is the tem- 
perature dependent contribution to the inverse square 
penetration depth). The temperature independent first 
term in the square brackets in the above equation, i.e. 
unity, should contain a factor S a p (from the first term 
in Eq. p4[ ). This factor S a p has been eliminated by im- 
posing the condition that the inverse square penetration 
depth should go to zero at T = T c where f{T c ) — 1. 

Eq. [35] for the penetration dfi»th can be seen to be 
identical to the result of LarkinEJ at T = 0, and to Eq. 
70 of LeggetttS at general T in the superconducting state. 
The model analyzed in this section is perhaps not, how- 
ever, identical to the models used in those articles. For 
example Leggett states that he considers a model that is 
isotropic, but at the same time states that translational 
invariance is not assumed while recognizing that this can 
only be an approximation. Thus, the basic assumptions 
of our work and those of Refs. [l?] and |l8|, while related, 
may be somewhat different. Fortunately the end results 
agree. 
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FIG. 1. The labelling of the nodes on the Fermi surface 
of YBa 2 Cu 3 06+^ 



VI. APPLICATION TO D-WAVE 
SUPERCONDUCTORS 

Now consider a cZ-wave high-temperature superconduc- 
tor modeled by a single copper-oxide plane having a 
Fermi surface such as that shown in Fig. |l|. From Eq. 

, which determines the Antisymmetric corrections to 
the quasiparticle energies, it is clear that only the values 
of 8e^ and at the Fermi surface nodes are relevant 
to the low energy properties (i.e. when ksT is much less 
than the maximum gap) . Also, the solutions of Eq. ^9| can 
be classified according to the irreducible representation 
of the point group C\ v (or 4mm) describing a tetrago- 
nal copper-oxide plane of a high T c superconductor, the 
independent solutions being 

Sei g = (Set + Sef + 8e£ + Sef) /4 
Sei y = (Sef - Set + 8e£ ~ Se£) /A 
8ei x = {5ef-5et-8ef + 5ei) /4 
8e£ y =(8ef + 8ef-8ef-8e£)/4 (36) 

where the indices 1,2,3 and 4 refer to the four nodes in the 
excitation spectrum, as defined in Fig. |l|. The external 
fields h£ at the nodes can be similarly classified. 
The solution of Eq. ^9] now gives the renormalized 

as 

h£ + 8e£{T) = h£/[\+Ff(T)] (37) 

with F^{T) = f^ln{2)k B T/(2irh 2 v F v 2 ). For d-w&ve su- 
perconductors, the quantity vp is the magnitude of the 
quasiparticle velocity at a node, as discussed following 
Eq. ^7] Also, r represents any of the irreducible repre- 
sentations present in Eqs. |3^. The /f^'s are defined by 

ft ^ A a i + + 2 A a 2 

rA. fa , ra c\ ra 

Jxy — 111 "T Jl3 — z /l2 

fi = fix - /iV (38) 

where /f 2 for example is f% k , for k and k' at nodes 1 
and 2, respectively. Note from Eq. [37] that a condition 
for the stability of the superconducting Fermi liquid is 
that [1 + Fjt(T)] > 0. Thus there is the possibility that 
the superconducting Fermi liquid will become unstable 
as the temperature is raised; this is perhaps of interest 
because of the non Fermi-liquid behavior of the normal 
state of the high T c superconductors. 

It is also useful to use Eq. to obtain an idea of 
how the (Symmetrical external fields are renormalized by 
Fermi-liquid interactions. First assume that there are 
no (Symmetrical external fields other than temperature, 
i.e. /if = (as is the case for the external magnetic 
fields of most interest in this article, which are purely 
Antisymmetrical) . Then the only term driving a nonzero 
contribution to Sef is the term on the right hand side 
proportional to f(E$ ). This term is proportional to 
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T 3 , thus giving a Sef oc T 3 , and will not be important in 
contributing to the properties of interest at the temper- 
atures satisfying fc^T <C A (A is the maximum gap). 
Thus the temperature-dependent contribution to Ssf will 

be neglected by putting f(E^) on the right hand side 
of Eq. ^9] equal to zero. In so far as the contributions to 
8ef driven by hf are concerned, it is clear from Eq. |28| 
that a knowledge of the Fermi-liquid interaction on the 
entire Fermi surface is required to calculate them, and 
that Sef must be determined on the entire Fermi surface 
(cf. Ref. ^l|). To obtain a rough idea of the nature of the 
solutions, consider a circular Fermi surface of radius kp 
and look for a solution of A g symmetry by considering 

a Fermi liquid interaction f^ k/ = f% , independent of k 
and k' , and a Symmetrical external field h 1 ^ indepen- 
dent of k. The solution, which is also independent of A; 
on the Fermi surface, gives the renormalized h A as 



h A +Se s A =h s A /[l + F 

f± a -Tin -Tin ' L 



(39) 



where F A 



-- kp j {i^fivp ). In contrast to the 

_4ntisymmetrical Fermi liquid parameters F^(T) ob- 
tained above, which go to zero linearly with temper- 
ature in the superconducting state in the clean limit 
(and hence have a dependence on temperature T ex- 
plicitly indicated) , the (Symmetrical Fermi liquid param- 
eter Fj[ is temperature independent and of approxi- 
mately the same magnitude as the corresponding nor- 
mal state Fermi liquid parameter. The same can be 
seen to be true of the Symmetrical Fermi liquid param- 
eters corresponding to other irreducible representations 
of Ci V . Note that the ratio of the _4ntisymmctrical to 
the Symmetrical Fermi- liquid F parameters is F A /F S » 
(f A /f s )[k B T/(hk F v 2 )}. 



A. penetration depth 

As noted above, the presence of a superfluid momen- 
tum contributes an _4ntisymmetrical external field to the 
Hamiltonian of Eq. This external field corresponds 
to the E irreducible representation of C± v with the p sx 
and p sy components of d, corresponding to the compo- 
nents Ex and Ey of Eq. |36|. The temperature-dependent 
quasiparticle contribution to the current density is thus 
easily evaluated using Eq. H with Eqs. |7|, |7] and ||. 



The result is 3 qp = r] qp p s where 



r)q P {T) = - 



2ln2e(v F ) 2 k B T 
[l+F E x (T)}irh 2 dv F v 2 



(40) 



Here, v F is the magnitude of the electrical transport ve- 
locity at a node, and d is the average spacing of the 
copper-oxide planes. Note that the Fermi liquid correc- 
tion does not alter the clean-limit linear-in-T contribu- 
tion to the r] qp (T), but rather makes a T 2 contribution 



(using (1 + F)^ 1 w (1 — F + ...)). Thus there are no 
external field renormalization corrections to the exper- 
imentally measured linear in T contribution to inverse 
square penetration depth. Hence the low temperature 
London penetration depth A is given by 



1 



1 



Al(T) A2(0) 



8ln2e 2 2 v F m 
c 2 h d v 2 



(41) 



where a = (v F /vp). 

For an isotropic Fermi surface and isotropic Fermi- 
liquid interactions, = v&(l + s-Ff), as described in 
the above discussion of isotropic s-wave superconductors 
(except for the factor | which arises due to our consider- 
ation here of a two-dimensional plane). By making use of 
this result, the result of Eq. [H], specialized to the case of 
an isotropic Fermi surface and Fermi- liquid interactions, 
can be seen to agree with the results of Refs. ^l] and |2^. 
A formula for the penetration depth for the anisotropic 
case has been derived independently in Ref. ^4|. In their 
result, our parameter v F is replaced by an expression de- 
pending on vp and on the Fermi- liquid interaction para- 
menters f%£, . The methods of this article can not make 
a connection between these two different expressions. 



B. Spin susceptibility 

The renormalization of the spin susceptibility due to 
Fermi-liquid interactions can be calculated in a similar 
way. The Zeeman interaction of the spin of an electron 
with the magnetic field contributes an Antisymmetric ex- 
ternal field of A g symmetry to the Hamiltonian. It fol- 
lows that the magnetic moment per unit area of a copper 
oxide plane is 



M 



Mb 
' L 2 



^2[f(E kA )-f(E ki _ 1 )]= X H 



where 



X(T) = 



Xo(T) 



l + FHry 



Xo(T) = 



fj? B ln2k B T 
7m vp v 2 



(42) 



(43) 



Note that here also the low-temperature clean-limit lin- 
ear in T magnetic susceptibility is not changed by Fermi- 
liquid interactions. These affect only terms of order T 2 
and higher in the susceptibility. 



C. Mixed-state density of states and specific heat 

In the clean limit for a d-wave superconductor, the 
quasiparticle density of states at zero energy in a mag- 
netic fieldEj of magnitude H normal to the copper-oxide 
planes varies as H 1 / 2 . Electron-electron interactions af- 
fect the magnitude of this H 1 / 2 contribution to the den- 
sity of states (and hence to the low-temperature specific 
heat) and this effect will now be calculated. 
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It follows from Eqs. [2?] and [37] that the zero-energy 
density of states in a single copper-oxide plane of a d- 
wave superconductor in a magnetic field (characterized 
by the superfluid momentum p s ) is 



N(0) = ( 2 



d 2 k 
(2n) 2 



<5(4 0) + (v* .p s )/[l + F#(T)]) 



4(K- Ps |) 



7rS 2 Wi,w 2 [l + F^(T)] 



(44) 



where is the value of \\ for k at a node. Here it is 
assumed that the superconductor is in the mixed state, 
and the angular brackets indicate an average of the super- 
fluid momentum p s over a unit cell of the vortex lattice. 
This density of states will give rise to a low-temperature 
specific heat varying linearly with temperature T. As 
in the evaluation of the low-temperature contribution to 
the penetration depth above, the quantity F F (T) (which 
represents the renormalization of the external field • p s 
by the Fermi- liquid interactions), can be neglected as it 
varies linearly with T and will contribute only to the 
next order T 2 contribution to the specific heai—pEq. 
can thus be seen to differ from previous workcZrEa only 
by the factor a = V f /vf, which takes into account the 
different effects of the electron-electron interaction on the 
electrical transport velocity and the quasiparticle veloc- 
ity. (This factor a is the same a that appears in the 
penetration depth formula, Eq. |il], above.) Thus, the 
equation for the clean-limit specific heat per unit volume 
in a magnetic field used in Rcf. [ll] (see also Refs. ^7] and 
E28|) becomes 



(45) 



The factor K is a numerical factor of order unity that 
depends on the way that p s (which is a function of po- 
sition in the .vortex lattice) is averaged over its position 
coordinate.cllES 

In the absence of a magnetic field, there are no Fermi- 
liquid corrections to the low-temperatiiee specific heat, 
which is given by the well-known rcsulttil 

C et = l8((3)k 3 B T 2 /(ndh 2 v F v 2 ). (46) 



VII. DISCUSSION AND CONCLUSIONS 

This article has developed an elementary phenomeno- 
logical approach to the evaluation of equilibrium Fermi 
liquid effects in superconductors with anisotropic Fermi 
surfaces and energy gaps. The approach is somewhat 
different from that used in a number of other articles, 
but gives a formula for the penetration depth that is in 
agreement with the results of Refs. [It], |lj, 2l] and ^ for 
the case of isotropic Fermi surfaces and Fermi-liquid in- 
teractions. This article extends these results to systems 




with anisotropic Fermi surfaces and Fermi-liquid interac- 
tions. A classification of the external fields in a supercon- 
ductor into what have been called above (Symmetric and 
Antisymmetric external fields was introduced. In <i-wave 
superconductors the (Symmetric external fields (including 
temperature) exhibit strong Fermi-liquid renormalization 
effects to which quasiparticles from the entire Fermi sur- 
face contribute (cf. Ref . |l|) , while the the Antisymmetric 
external fields (including the magnetic fields that act on 
both the orbital and spin degrees of freedom of the quasi- 
particles) exhibit relatively weak temperature-dependent 
renormalizations that arise from the nodal quasiparticles 
only, and that can often be neglected. 

Eq. 41 above for the London penetration depth has 
exactly the same form as Eq. 6 of Ref. [yj which finds 
(from a detailed analysis of a number of experiments) 
a 2 = (v F /v F ) 2 = 0.43 for Bi 2 Sr 2 CaCu 2 8 and a 2 = 
{v f /vf) 2 = 0.46 for YBa 2 Cu30j_5. This implies a value 
of the electrical transport velocity v F at a nodal point 
significantly smaller than quasiparticle velocity vf- 

To appreciate the significance of this result concerning 
the ratio of v F to vp, some idea of the expected relative 
magnitudes of v F and v F must be obtained. Recall that 
the model studied in this article is that of a single band 
of electrons interacting with each other and with a pe- 
riodic potential. All of the interactions of the electrons 
in this single band with electrons in other bands or with 
the ion cores are incorporated into the underlying peri- 
odic potential. The solution of the problem of an electron 
in this periodic potential in the absence of any electron- 
electron interactions gives what was called the bare band 
energy s\ and the bare band velocity w\ = Wet / dk. At 
this level of approximation the quasiparticle velocity and 
the electrical transport velocity are the same, and equal 
to the bare band velocity. Now if electron-electron in- 
teractions are added in the Hartree-Fock approximation, 
the band energy, the Fermi surface, and the quasiparticle 
velocity will all be modified, but the electrical transport 
velocity remains equal to the bare band velocity. Further- 
more, if the electron-electron interactions are considered 
in a higher order of approximation in an isotropic model, 
there will be further changes to the quasiparticle velocity, 
but the electrical transpoct,vclocity will still remain equal 
to the bare band velocityEj. Two studies of the effective 
mass ratio m*/m in an electron-gas model in a high- 
density approximation have been summarized in Fig. 10 
of Ref. |2{| (For an isotropic electon gas model the effec- 
tive masses are defined by Vp = pf/tti* and v F — p F /m.) 
These studies show that for high densities to* /to is less 
than unity by up to 5%, but when the density is decreased 
the ratio of m* /to becomes greater than unity, by up to 
10% for the range of parameters studied. Thus, at lower 
densities, where stronger correlation effects are expected 
(because the ratio of the potential energy to the kinetic 
energy in the electron gas becomes larger), the quasipar- 
ticle becomes heavier, i.e. tends to move more slowly, 
as would naively be expected. It is not easy to fit the 
experimental results on high T c superconductors quoted 
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above into the framework of these results. The experi- 
mental result implies an effective mass ratio m* jm less 
than unity. Effective mass ratios less than unity (but 
by a maximum of 5%) are found in Ref. ^9], but only in 
the higher density and relatively less strongly correlated 
electron gases. This does not seem to fit very well with 
the idea of electron correlations that are strong enough 
to lead to an antiferromagnetic Mott insulating state in 
the high T c materials for certain doping levels. It is prob- 
ably to be expected that high-density electron gas ideas 
are not particularly useful in understanding the high T c 
materials. 

Another potentially relevant study is that of Ref. |3(], 
which has calculated the ratio of the quasiparticle ef- 
fective mass to a bare band mass in a model relevant 
to high T c superconductors (see their Eq. (9) and sur- 
rounding discussion). However the band mass there is 
envisaged as being derived from a energy band struc- 
ture closely related to that determined by photoemission 
spectroscopy, and presumably already includes some ef- 
fects of the electron-electron interactions. Thus it is not 
clear how this quantity might be related to the electrical 
transport velocity which is required for an interpretation 
of the penetration depth measurements in the formula 
given above. It is nevertheless of interest to note in Fig. 



5 of Ref. 30 



that electron interaction effects produce a 
slight band broadening, which goes in the right direction 
of accounting for the observed value of a. 

Fermi-liquid corrections to the spin susceptibility, and 
to the magnetic field induced zero-energy density of 
states have also been evaluated above for d-wave super- 
conductors at temperatures much less that T c . 
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